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Abstract 

We have developed a two-dimensional model of quantised vortices in helium II moving under the 
influence of applied normal fluid and superfluid in a counterfiow channel. We predict superfluid 
and vortex-line density profiles which could be experimentally tested using recently developed 
visualization techniques. 

PACS numbers: 67.25.dk, 47.37.+q 
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I. MOTIVATION 

The recent development of visualization techniques in superfluid helium based on micron- 
size tracers has raised the possibility of experimentally determining superfluid and 
normal fluid profiles and the spatial distribution of quantised vortices in a channel, thus 
solving outstanding problems in quantum turbulence. Particularly interesting (and relevant 
to engineering applications) is the turbulence induced by heat transfer (counterflow turbu- 
lence). It is well known that if the applied heat flux Q is less than a small critical value Qc 
then the heat is carried by the normal fluid component, and the superfluid component flows 
in the opposite direction to conserve mass, li Q > Qc the superfluid component becomes 
turbulent, forming a disorganized tangle of quantised vortices. In some geometries, at larger 
heat flux a transition to a more intense vortex tangle has been observed |3|, which is perhaps 
related to the onset of turbulence in the normal fluid ^. 

In this paper we are concerned with the intermediate heat transfer regime, in which the 
normal fluid is still laminar, but the superfluid forms a turbulent tangle. The issue which 
we address is the average superfluid profile and the average spatial distribution of vortices 
in the channel. Let v„ and be respectively the normal fluid and superfluid velocity flelds. 
The normal fluid satisfles no slip boundary conditions at the wall of the channels, so, for the 
sake of simplicity, we assume that v„ is a steady classical parabolic Poiseuille proflle (this 
approximation implies that the constant pressure drop along the channel is the same with 
and without the vortex tangle). The superfluid slips at the boundaries, so, in the absence of 
vortices, it is natural to assume that has a constant (uniform) proflle. The velocity flelds 

and are related by the counterflow condition of no net mass flow along the channel. 

The question which we ask is the following: if Q > Qc, are the vortices uniformly dis- 
tributed in the channel, or do they organize themselves spatially, creating a non-uniform 
superfluid proflle at scales larger than the average vortex separation i but smaller than the 
channel size D ? 



II. MODEL 

To answer the question we consider the following idealised two-dimensional model which, 
we argue, captures the most important physical ingredients. Let x and y be respectively 
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directions along and across the channel, with walls at y = ±D/2 and periodic bound- 
ary conditions at x = and x = X. The normal fluid velocity is v„ = [v^, v^) = 
(— Kio[l — (2y/-D)^] 7 0) with Vno > and thus pointing in the negative x direction. The 
superfluid velocity = {v^, v^) is decomposed in two parts, v., = v^q + v^j; the former 
is the uniform flow v^o = Ko(^)x , x being the unit vector along x, and the latter, v^j, is 
the velocity field induced by vortex points located at positions = {xj{t),yj(t)) (for 
j = 1, ■ ■ ■ A^) where t is time. Half the vortices have positive circulation Tj = k, and half 
have negative circulation Tj = —k, where k, = 10^^ cm^/s is the quantum of circulation in 
superfluid ^He. To satisfy the superfiuid's boundary condition that = at y = ±D/2 (no 
flow into the wall), we attach to each vortex point an infinite series of image vortex points 

n 

in the positive and negative regions y > D/2 and y < —D/2 [5j. 
The equation of motion of a vortex located at rj is j6| 

dv ■ 

= Vs(rj) + as'{rj) x (v„(rj) - v,(rj)) + a'(v„,(rj) - v,(rj)), (1) 

where s^- is the unit vector along the vortex j (in the positive or negative z direction) and a 
and a' are temperature dependent mutual friction coefficients 0]. 

The quantity Ko is determined at each time t by imposing the counterflow condition 

p„ < < > +ps < < >= 0, (2) 

where p„ and ps are the temperature-dependent normal fluid and superfluid densities, p = 
pn + Ps is the total helium density, and < > and < > are the channel averages of the 
x-components of the normal fluid and superfluid velocities, defined as 



1 /-A fD/2 I rX rD/2 

<^s>=Y7^ / v^ix,y)dxdy = Vso + / v'li{x,y)dxdy (3) 

XDJq J-d/2 XDJo J~d/2 



and 

1 /■^/2 2 
< < >= n / , <Wy = -^Ko- (4) 

D J~D/2 6 

To model the creation and the destruction of vortices within our 2-dimensional model, we 
proceed as follows. When the distance between two vortex points of opposite circulation 
becomes smaller than a critical value ei, we perform a "numerical vortex reconnection" 
and remove these vortex points; similarly, when the distance between a vortex point and 
a boundary is less than €2 = 0.5 ei, we remove this vortex point (the vortex of opposite 
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circulation being the nearest image vortex beyond the wall). To maintain a steady state, 
when a vortex point is removed, a new vortex point of the same circulation is re-inserted 
into the channel, either on the axis [y = 0) or near the walls or randomly, as described in 
the next section. 

At selected times, given the vortex configuration rj(t) (j = 1, ■ ■ ■ A^), we define a coarse- 
grained superfluid velocity by averaging the components of the (microscopic) velocity 
Vs over channel strips of size A in the y direction, such that i < A < D. The limit 
i << A << D corresponds to the Hall - Vinen - Bekarevich - Khalatnikov (HVBK) equations 
jsj. We assume that the velocity is parallel to the walls, = {vg, 0). The curl of can 
be interpreted as the coarse-grained superfiuid vorticity. 

To make connection with the experiments, we interpret n = N/ (XD) (number of vortex 
points per unit area) as the vortex line density L (vortex length per unit volume), from which 
i = n~^/^ is the average intervortex spacing. The imposed heat fiux Q, which is reported in 
experiments, is related to the average counterflow velocity Vns =< > — < > hy the 
relation Q = TpsSVns, where T is the temperature and 5* the specific entropy. 

The calculation consists in computing the evolution of the vortex points, starting from 
an arbitrary initial condition, until a steady state regime is achieved and the profile of 
becomes constant. The time integration is performed using the second-order Adam-Bashfort 
method with time step At. 

III. RESULTS 

The numerical code solves the governing equations written in dimensionless form. 
The units of length, velocity and time are respectively 6c = D/2 = 4.55 x 10~^ cm , 
Uc = K,/{27i6c) = 3.49 X 10~^ cm/s, tc = Sc/uc = 0.13 s. Non-dimensional quantities are 
denoted by the superscript We choose parameters taking into account the available com- 
puting power and the experiments of Tough and collaborators N = 1876, A* = 6, 
T = 1.7 K, D* = 2 (corresponding to tube i?4 in Ref. [sl), ei = 1.25 x IQ-^^, A = 1.25 
V*g = —478.5 and At* = 1.9 x 10~^. At this temperature, taking into account pressure and 
temperature variations experimentally measured along the channels j^, we can assume the 



constant values = 0.112 g cm^^ p„ = 3.32 x 10"^ g cm^^ and S = 0.395 J (g K)-^ (ll|. 
Using Eqs. ([2]), dlD and the definition of Vns "we have V*q = 553.6. The vortex density is 
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n* = 156.3, and the average intervortex spacing, i* ^ 0.08, corresponds to the dimensionless 
number L^/'^D = D*y/n* = 25 , which is typical of counterflow experiments 10|. 




FIG. 1: (color online). Initial random vortex configuration (left) and vortex configuration at 
t* = 4.75 X 10^^ (right). Positive and negative vortices are respectively denoted by empty (red) 
and filled (black) circles. 

To distinguish the three different "re-nucleation" regimes, hereafter we indicate with (a) 
random vortex-reinsertion in the channel, and with (b) and (c) vortex re-insertion on the 
axis and near the walls, respectively. 

Case (a) is the simplest and corresponds physically to 3-dimensional superfluid vorticity 
production taking place throughout the vortex-tangle due to reconnections and vortex-rings 
emission. The initial condition, shared by all three cases, is shown in Fig. [T] (left) and 
consists of a random vortex configuration. During the evolution of case (a), positive and 
negative vortices move on the average towards the y* = —1 and y* = +1 wall respectively. 
The trajectory of an individual vortex can be very irregular, due to the interaction with 
other vortices. A drift in the positive x direction is superimposed to this motion towards 
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walls. 

After a time- interval of the order of r* ^ 2 x 10~^, the vortex configuration reaches a 
steady state, as illustrated in Fig. [T] (right). Note that there are more positive vortices in 
the y* < region and more negative vortices in the y* > region, as confirmed by the 
density profile of positive vortices, n^{y*), shown in Fig. [3] (right) (the density profile of 
negative vortices, n*_{y*), is symmetrical with respect to the channel axis). The polarization 
of the vortex configuration is not complete - not all positive vortices are in < and not 
a// negative vortices are in y* > - in agreement with arguments discussed in Ref. 12 1. 

The resulting coarse-grained superfiuid velocity profile associated to this partial polar- 
ization is compared to the driving normal fluid velocity in Fig. H] (left): we find that the 

2 19 , 

profile is almost parabolic, v* y* ' . The distribution of vortices irrespective of their sign, 
n*{y*), is almost constant, as shown in Fig. |3] (left). 

The other two cases which we investigate, case (b) (nucleation at y* = 0) and case (c) 
(nucleation at y* = ±1), are suggested by the analysis of Geurst [l3| based on the HVBK 
equations. Case (b) induces a vortex density profile n*{y*) with a sharp peak near y* = 
(Fig. |3] (left) ). We find that this concentration of vortices on the axis is unstable: collective 
motion of vortices of the same polarity appears, and the resulting coarse-grained superfiuid 
velocity profile, v*{y*), shown in Fig. |2] (right), undergoes small-amplitude but persistent 
oscillations. A power-law regression yields the dependence v* ~ y*^ ^] as in case (a), after 
a transient of the order of r^* ~ 6 x 10"^, we obtain the partially polarized steady state 
vortex-distribution shown in Fig. |3] (right). 

Finally, we consider case (c). This nucleation regime corresponds physically to 3- 
dimensional superfiuid vorticity production arising from vortex-rings emission by vortices 
pinned to the walls. In our model, the positive (negative) vortex points are re-inserted near 
(rather than at) the y* = 1 {y* = —1) wall, in order not to collapse rapidly on the boundaries 
due to the Magnus-mutual friction forces balance. Let ^ be the distance of nucleation away 
from the walls. To investigate the dependence of the flow on ^, we choose values C,i = i and 
^2 = 0.5 £, which we refer to as cases (ci) and (02). The resulting coarse-grained superfiuid 
velocity profile v*{y*) in the two cases is shown in Fig. |2] (left): note the reversed concavity, 
compared to cases (a), (b) and the normal fluid. 

Fig. [3] shows that the nucleation location clearly influences the vortex distribution and 
the steady state vortex density profile n*{y*), which is reached after a transient time interval 
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FIG. 2: (color online). Solid (red) line: coarse-grained superfluid velocity v*{y*); the labels (a), 
(ci), (C2) (left) and (6) (right) correspond to the cases discussed in the text. Dashed (blue 
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line: analytical laminar solution [1^ of HBVK equations which applies Geurst's approach 
to Cartesian geometry. Curves Gc (left) and Gb (right) correspond respectively to cases (c) and 
(b) in the text. Dot-dashed (black) line: normal fluid velocity profile (left, right). 



of the order of r* 3 x 10 ^, independently of ^. 



IV. DISCUSSION 



The model which we have presented, being 2-dimensional, is clearly rather idealized, and 
we do not claim that it is possible to make direct comparison with experiments. Nevertheless, 
the model contains what we argue are the most important physical ingredients, and it allows 
us to make predictions about the vortex distribution and the superfluid profile across a 
channel at small values of the applied heat flux. The 2-dimensional solution which we have 
found can be interpreted in 3-dimensions as vortex loops which move from the outer parts of 
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FIG. 3: (color online). Vortex density n*{y*) (left) and positive vortex density n^(y*) (right) for 
cases: (a), dashed (blue) line; (b), dot-dashed (black) line; (ci) and (02), solid (red) lines. 



the channel towards the centre, speeding up during this process, or, perhaps more precisely, 
as a tangle which is polarised by the presence of such loops. 

Our coarsed-grained superfluid velocity profile v*{y*) compares very well with the 2- 
dimensional laminar solution of the HVBK equations which can be derived by applying 
Geurst's approach 13| for a cylindrical pipe to Cartesian geometry [l^. Indeed, Fig. [2] 
(right) shows that case (b) and this laminar solution jl^ are very similar (with the exception 
of the channel axis, where Guerst's solution is singular and n* — 00). 

In Fig. [2] (left) we compare the analytical solution corresponding to nucleation at the 
walls {C, = 0) to the velocity profile v*{y*) for case (c). It is apparent that, with decreasing 
^, Vgiy*) tends to the analytical solution (the difference which is present near the walls arises 
from the boundary conditions for the superfluid velocity which ensure infinite vorticity in 
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the nucleation region 

,Q)- 

Further work will generalise the approach which we have presented, by solving self- 
consistently the equation for the normal fluid in the presence of the mutual friction, rather 
than assuming a given profile for v^. We also plan to investigate the dependence of the 
profiles on the vortex densities. The current rapid progress of visualization techniques, such 



as the recently developed laser-induced fluorescence of metastable molecules 15|, will clearly 



stimulate more work on the nature of laminar and turbulent profiles of helium II in channels. 
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